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Abstract 

The complete next-to-next-to-leading order (i.e. 0(v 2 ), 0(va s ) and 0(a 2 )) relativistic cor- 
rections to the total photon mediated tt production cross section at threshold are presented 
in the framework of nonrelativistic quantum chromodynamics. The results are obtained 
using semi-analytic methods and "direct matching" . The size of the next-to-next-to-leading 
order relativistic corrections is found to be comparable to the size of the next-to-leading 
order ones. 



PACS numbers: 14.65.Ha, 13.85.Lg, 12.38.Bx. 



Due to the large top quark mass, which allows for the decay channel t — ► bW, ti production at 
threshold in lepton pair collisions offers the unique opportunity to study heavy quark-antiquark bound 
state and near threshold dynamics using perturbative QCD Q. With this motivation in mind a 
considerable number of studies has been carried out in the past in order to calculate tt production 
observables |2|, |3|, ||, ||, || and explore their potential for measurements of the top quark massQ M% and 
the strong coupling a s at future experiments like the LC (Linear Collider) Q or the FMC (First Muon 
Collider) ||]. In view of the high precision which might be achieved for the QCD calculations as well as 
for ti production measurements even relatively small effects coming from a light Higgs boson g [9| have 
been investigated. However, the present day analyses only include QCD effects up to next-to-leading 
order (NLO) in form of the one-loop corrections to the QCD potential (TO, 11] and various 0(a s ) short- 



distance corrections. A complete next-to-next-to-leading order (NNLO) calculation, which would be 
necessary to study the reliability of the present day analyses and to make the consideration of small 
effects from beyond QCD at all feasible, has been missing so far. 

In this letter we present the complete NNLO relativistic corrections to the total photon medi- 
ated tt production cross section. As NNLO we count all corrections of order v 2 , va s and a 2 relative 
to the cross section in the nonrelativistic limit, v being the cm. velocity of the top quarks. As rela- 
tivistic, on the other hand, we count the corrections coming from the top quark kinetic energy, the 
ti production and annihilation process including the short-distance corrections, and the tt interaction 
potentials. We use NRQCD [12] to conveniently parameterize calculations and results in a systematic 



manner following the approach proposed in [O]. The calculations are carried out using semi-analytic 



methods and the "direct matching" procedure introduced in [14]. We would like to point out that 
NNLO corrections involving the top quark decay are not determined here. The latter effects would 
include 0(a 2 ) two loop corrections to the free top quark width and a consistent treatment of the ef- 
fects from the off-shellness of the decaying top quarks, the time dilatation and the interactions among 
the decay products and the other top quark (if it is not decayed yet). Although the size and the 
interplay of all these effects have been studied at various places in the literature (see e.g. || [l5|, |i~6f), 
their consistent treatment at NNLO still remains an open problem. As far as the NNLO relativistic 
corrections discussed in this letter are concerned we will use the naive replacement 

E = y/s-2M t — >E = E + iT t (1) 

in the spirit of |IJ in order to examine their size and properties, where Tt represents a constant which 
is not necessarily the decay width of a free top quark. We also would like to emphasize that we 
treat all interactions purely perturbatively and that nowhere in this work the confining long-range 
contributions to the QCD potential or other nonperturbative effects are taken into account. This is 
somewhat contrary to the standard present day approach used to describe tt production at threshold 
(see g, H ||), but we take the position that nonperturbative effects might be added later as a 
correction. 

We start from the NRQCD Lagrangian 
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1 Throughout this paper Mt is understood as the top quark pole mass. 
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The gluonic and light quark degrees of freedom are described by the conventional relativistic La- 
grangian, whereas the top and antitop quark are described by the Pauli spinors ijj and Xi respectively. 
For convenience all color indices are suppressed. The straightforward antitop bilinears are omitted 
and only those terms relevant for the NNLO cross section are displayed. Dt and D are the time and 



space components of the gauge covariant derivative -D», and E l = G° l and B l 



and magnetic components of the gluon field strength tensor (in Coulomb gauge). The short-distance 
coefficients ay, . . . , 05 are normalized to one at the Born level. Because we use "direct matching" [14] 
their higher order contributions are irrelevant for this work. 

To formulate the normalized total tt production cross section (via a virtual photon) R = 
<r(^ + ^_ — > 7* — > tt)/a p t ((Jpt = 4-7ra 2 /3s) in the nonrelativistic region at NNLO in NRQCD we start 
from the fully covariant expression for the cross section 



R(q 2 ) = ^-Im[(0\TUq)~f(- 



0: 



(3) 



where Qt = 2/3 is the top quark electric charge. We then expand the electromagnetic current (in 
momentum space) jn(±q) = (^y^t) (±5) which produces/annihilates a tt pair with cm. energy \f~q^ in 
terms of 3 Sy NRQCD currents up to dimension eight (i = 1, 2, 3) 



l(q) = h (q) ~ ^ (^(-§ D)x) (q) + 



(4) 



where the constants b\ and hi are short-distance coefficients normalized to one at the Born level. The 
expansion of ji(—q) is obtained from Eq. (0) via charge conjugation symmetry. It should be noted 
that only the spatial components of the currents contribute. Inserting expansion (||) back into Eq. (||) 
leads to the NRQCD expression of the nonrelativistic cross section at the NNLO level 
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where 
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To obtain the factor ^irQl / 3M t 4 in the second line of expression (||) we have already used relation (|l4| ) . 
The cross section is expanded in terms of a sum of absorptive parts of nonrelativistic current correla- 
tors (containing long-distance physics^) multiplied by short-distance coefficients Cj (i = 1,2,.. .). In 
Eq. (|H|) we have also shown the dependences on the various renormalization scales: the soft scale fj, so &, 
and the hard scale //hard are governing the perturbative expansions of the correlators and the short- 
distance coefficients, respectively, and arise from the light degrees of freedom in £ N rqcd,0 whereas 
the factorization scale /if ac essentially represents the boundary between hard (i.e. of order Mf) and 

2 In the context of this paper "long-distance" is not equivalent to "nonperturbative" . 



Throughout this work we use the convention a s — a 



("(=5) 



in the MS scheme. 
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soft momenta. This boundary is not uniquely denned and therefore both, the correlators and the 
short-distance coefficients, in general depend on /if ac (leading to new anomalous dimensions). Because 
the term in the second line in Eq. (g) is already of NNLO (i.e. suppressed by v 2 ) we can set C2 = 1 
and ignore the factorization scale dependence of the correlator Ai- The calculation of all terms in 
expression (|5|) proceeds in two basic steps. 

Step 1: Calculation of the nonrelativistic correlators. - The calculation of the correlators at NNLO is 
simplified enormously by the fact that for the ti pair, which is produced and annihilated in a color sin- 
glet S-wave state, all retardation effects can be neglected. Technically this means that the transverse 
gluon exchange, which would lead to temporally retarded interactions (closely related to Lamb-shift 
type corrections known in QED), can be treated as instantaneous, i.e. its energy dependence can be 
ignored. This can be seen by using formal counting rules (see e.g. [17, or from physical con- 
siderations because a real gluon (i.e. one that carries energy) which is radiated after the ti pair was 
produced has to be absorbed before the ti pair is annihilated. Phase space effects and the dipole 
matrix element then lead to a suppression of this process oc t> 3 , i.e. to an effect beyond NNLO. In 
other words, as far as the nonrelativistic correlators in Eq. (||) are concerned, NRQCD reduces to a 
two-body (top-antitop) Schrodinger theory. The potentials in the resulting Schrddinger equation are 
determined by considering ti — > ti one gluon exchange t-channel scattering amplitudes in NRQCD. To 
NNLO (i.e. including potentials suppressed by at most a 2 , a s /Mt or 1/M 2 relative to the Coulomb 
potential) the relevant potentials read (a s = a s (fi so { t ), C A = 3, Cf = 4/3, T = 1/2, jl 
r = \r\) 
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where St and Sj are the top and antitop spin operators and L is the angular momentum operator and 
(ni = 5) 
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The constants (5q and (5\ are the one- and two-loop coefficients of the QCD beta function and 
7 = 0.577216 ... is the Euler constant. V c is the Coulomb (static) potential. Its 0(a s ) and 0(a 2 s ) cor- 
rections have been determined in []To| , p"l| 1 and [jl9| , respectively. Vbf is the Breit-Fermi potential known 
from positronium and V NA a purely non-Abelian potential generated through non-analytic terms in 
one-loop NRQCD (or QCD) diagrams containing the triple gluon vertex (see e.g. |20| | for older refer- 
ences). The nonrelativistic correlators are directly related to the Green function of the Schrodinger 
equation 



V 2 



V' 1 



+ V c (r) + V BF (r) + V NA (f) - E ) G(f, f, E) = 6® (f-r 



where Vbf is evaluated for the 3 Si configuration only. The correlators read 

Ai = 6N C 
A 2 = M t EAi. 
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(12) 

(13) 
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Relation (13) can be easily inferred by taking into account that the Green function G(f, r 1 , E) describes 
the propagation of a top-antitop pair which is produced and annihilated at distances \f\ and |f |, 
respectively. Because the exact solution of Eq. (12) seems to be an impossible task, we rely on a 
numerical solution of the equation 



+ V c (r)-E^jG c (r,f*,E) = ff( 3 )(f-^) 



(15) 



using techniques developed in 0, Q. The result for G c (0, 0, E) is then combined with the corrections 
to the leading order (LO) Coulomb Green function G^ [^] (defined through Eq. ( |T5| ) for V c (r) = 
—C F a s /r) coming from the kinetic energy correction and the potentials Vbf and Vn A . These corrections 
are calculated analytically using Ray leigh- Schrodinger time-independent perturbation theory, 



5G(f, r, E) 



dx 3 G^\r,x,E) 



V 4 

4M 4 3 
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(16) 



For the calculation of Eq. fllq ) we use techniques employed in p2| , where the Abelian NNLO contri- 
butions have already been determined. The final result for A\ at NNLO reads 
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and \& is the digamma function, ^f(z) = 4- lnT(z). In the first line of Eq. (|T7|) the LO Green function 
has been subtracted to avoid double counting of the LO contribution contained in the third line. It 
should be noted that the limit |r| , |f| — ► in expression ( |i~6l) causes UV divergences which are regulated 
using the short-distance cutoff ^f ac . Further, all power divergences oc [i{ ac /M t are subtracted |23|] and 
jUfac is defined in a way that expression ( |T7| ) takes the simple form shown above. The corresponding 
result in any other regularization scheme could be obtained from the one presented here through a 



redefinition of the factorization scale. For A2 only the LO contribution in Eq. (17) is relevant and we 
arrive at 



A 2 
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There are no non-Abelian contributions to A%- 

Step 2: Matching calculation. - The contributions to C\ up to 0(a 2 ) are determined by matching 
expression (||) directly to the two- loop cross section calculated in full QCD in the (formal) limit 
a s <C v <C 1 for stable quarks (1^ = 0) keeping terms up to a 2 and NNLO in the expansion in v and 
setting // so f t = /ihard-0 To obtain the contributions oc a 2 originating from the potential V c in this limit 
we also employ time- independent perturbation theory in analogy to Eq. (|l6|). The corresponding cross 
section in full QCD reads (ah = a s (/^hard)> v = (E /M t ) l l 2 ) 



pNNLO 
-"^loop QCD 



N C Q 



3 

- v 
2 



17 



+ 



C 2 F 7T 2 3^ 

— + 2 Cf 



+ 



A9C 2 F 7T 2 



16 

2Cf + Ca 



7T 



3vr^ 



7T 5 
6 V H 1> 

2 



192 



7T" 



3 C F /11_ 
+ K + ~4 hr CU 



11, 4w 2 M? 
— In — ^ — - + 

24 /^hard 

2Tn,) ln^pL 

A'hard 



31 
72 



, Av 2 M? 

6 Kard 



_5_ 
18 



(20) 



where 
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The Born and 0(a s ) [p4|| contributions are standard. At order a 2 the contributions in Eq. (20) 
proportional to C F , CaCf, CfTh\ and CfT have been calculated in |2 
respectively. The result for C\ reads 



and 21 2t 
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The consistency of the "direct matching" procedure ensures that C\ does not contain any energy- 
dependent terms. We would like to point out that the factorization scale dependence of Im[«4i] is 



cancelled by the factorization scale dependence in C\ up to a small term oc Cp ^siu" m (u^") ( see 



M t 



also [22]). This term remains as a consequence of our ignorance of a consistent treatment of the finite 



For a s <C"<l, i.e. far away from the threshold regime, no distinction between soft and hard scales is needed. 
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Figure 1: (a) The total normalized photon-mediated tt cross section at LO (dotted lines), NLO 
(dashed lines) and NNLO (solid lines) for the soft scales /i so ft = 50 (upper lines), 75 and 100 GeV 
(lower lines), (b) The NNLO cross section for a s (M z ) = 0.115 (solid line), 0.118 (dashed line) and 
0.121 (dotted line). More details and the other parameters are given in the text. 



width effects at the NNLO level. Due to the small size of this contribution, however, the corresponding 
ambiguity can be ignored for the examination of the relativistic NNLO corrections determined in this 
paper. The LO cross section can be recovered from the NNLO one in Eq. (^) by taking into account 
only the dominant contributions in the third line of Eq. (17) and setting C\ = 1, C2 = 0, whereas 
the NLO cross section can be obtained by incorporating also the 0{a s ) corrections to the Coulomb 
potential and to the constant C\. In Fig. |l|(a) the LO (dotted lines), NLO (dashed lines) and NNLO 
(solid lines) normalized cross sections are plotted versus E in the range —5 GeV < E < 5 GeV for 
M t = 175 GeV, a s {M z ) = 0.118 and T t = 1.43 GeV. For the scales the choices /x soft = 50 (upper 
lines), 75 and 100 GeV (lower lines) and /Uhard = Mfac = Mt have been made and two- loop running of 
the strong coupling has been used. It is evident that the NNLO corrections are large. Compared to 
the NLO cross section, the 15 peak is shifted towards smaller energies by several hundred MeV and 
the large negative NLO corrections for positive energies are compensated to some extent. Whereas 
the location of the 15 peak is quite insensitive to changes in the soft scale, the residual dependence of 
the normalization of the NNLO cross section on the soft scale /x so ft is not improved at all compared to 
the NLO cross section. In fact, for energies above the 15 peak it is worse for the NNLO cross section 
than for the NLO one. The dependence of the NNLO cross section on the hard scale //hard an d the 
factorization scale /Xf ac are, on the other hand, much smaller and, therefore, not displayed here. The 
behavior of the NNLO corrections clearly indicates that the convergence of the perturbative series 
for the tt cross section is much worse than expected from the general arguments given by Fadin and 
Khoze [0. For the normalization of the cross section we estimate, at least at the present stage, a 
theoretical uncertainty at the level of five to ten percent. For comparison, in Fig. |l|(b) the NNLO 
cross section is displayed for a s {M z ) = 0.115 (solid line), 0.118 (dashed line) and 0.121 (dotted line) 
and /i S oft = 75 GeV. The other parameters are chosen as before. A more detailed examination of the 
NNLO contributions will be carried out in a future publication. 

This work is supported in part by the U.S. Department of Energy under contract No. DOE DE- 
FG03-90ER40546. 
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